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Abstract 

A dilaton potential is adjusted to recently confirmed lattice QCD thermodynamics data in the 
temperature range (0.7... 3.5)T C where T c = 155 MeV is the pseudo-critical temperature. The 
employed holographic model is based on a gravity-single-field dilaton dual. We discuss conditions 
for enforcing (for the pure gluon plasma) or avoiding (for the QCD quark-gluon plasma) a first-order 
phase transition, but still keeping a softest point (minimum of sound velocity). 
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I. INTRODUCTION 


The celebrated AdS/CFT correspondence [T] has sparked a large number of dedicated in¬ 
vestigations of strongly coupled systems (cf. [2] for recent surveys). A particular field of 
application is provided by the strong coupling nature of QCD at low momentum/energy 
scales. While employing 5d Einstein gravity in the dual description is strictly justified only 
in the large N c and large t’Hooft coupling limits of the boundary theory, which ensure 
suppression of loop and stringy corrections to classical gravitation theory, by such means 
nevertheless one could study models which are expected to exhibit a behavior resemblant 
of QCD. The aim is then often to understand, on a qualitative level, phenomena which are 
hardly accessible in the 4d quantum field theory. A prominent example is given by real¬ 
time phenomena, e.g. within QCD. Other phenomena, such as the hadron spectrum or the 
equation of state, are accessible by lattice QCD calculations - but here one would like to 
understand qualitatively the emerging numerical results by means of transparent models. 
Many facets of the QCD equation of state are fairly known by now, both for the physical 
parameter section and for various limits of parameters (e.g. quark masses, dimension of 
the gauge group, flavor number, adjoint representations of quarks etc.). This statement 
applies only for finite-temperature (T) QCD at zero baryo-chemical potential (//). However, 
in relativistic heavy ion collisions, the bulk of excited matter has fi > 0, as inferred from 
the analysis of hadron abundancies [3]. The knowledge of the QCD equation of state is, 
therefore, presently incomplete (in particular beyond the range accessible by the /i/T <C 1 
expansion) and calls urgently for an improvement. In particular, there are several ideas that 
QCD allows for a critical point in the T — /j phase diagram where the cross-over turns in a 
first-order phase transition. Mainly based on universality arguments, a multitude of models 
have been employed to locate the critical point ®Ej, but also more directly QCD anchored 
approaches, e.g. Dyson-Schwinger equations as integral formulation of QCD, have been used 
[6]. Parallel to the theoretical attempts, also special experimental searches are conducted, 
e.g. the beam energy scan at RHIC [7]. 

Coming back to options for modeling a phase diagram similar to QCD with the conjectured 
critical point, we mention [8j, where in a holographic model, including gravity, a dilaton 
field and a U(l) gauge field, the possibility of such a realisation has been demonstrated. 
The set-up of PJ is based on a dilation potential, which features qualitatively the equation 
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of state at n — 0, supplemented by a dynamical strength function, which is adjusted to 
the quark number susceptibility, again at /i = 0 . While in the infinitely heavy quark mass 
(' m q —> oo) limit of QCD, which becomes then a pure Yang-Mills theory, the equation of 
state is known since some time [9] and has been confirmed by high-precision lattice QCD 
simulations [10], the status of QCD with physical quark masses has been settled only very 
recently. After refinements in the lattice discretization schemes and actions and continuum 
extrapolations the results of two independent collaborations HD El became consistent. 
Given this new situation and having in mind e.g. an application in the spirit of [5] to the QCD 
phase diagram modeling, one should seek for an appropriate dilaton potential, reproducing 
sufficiently accurately the by now known QCD equation of state at // = 0. This is the aim 
of the present note. We have hereby the attitude to take the AdS/CFT dictionary literally, 
i.e. translate, without corrections due to N c = 3 or finite coupling, the 5d Riemann metric 
into 4d energy-momentum tensor components (or correlators) and vice versa. 


II. ADJUSTING A DILATON POTENTIAL 


At /i = 0, the equation of state, in parametric form, follows from 


f <t>H 


LT((j> H ) = exp (A(0 0 ) + / d0 




G 5 s{<j> H ) = ^ ex P (3A(0o) + - J 


f <t>H 


for entropy density s and temperature T, where the scalar function X(0;0#) 
mined by the system (a prime means a derivative w.r.t. 0) 


( 1 ) 

( 2 ) 

is deter- 

( 3 ) 

( 4 ) 


which is integrated from the horizon <pn — e to the boundary 0o with initial conditions 

• Y( ^- e) = -ivHf| + 0(tl) - (5) 

Y(4h - c) = -- C MlLzA + O(e 0 ), ( 6 ) 

and e —> 0. The quantity A(4>o) encodes the near-boundary behavior of the model. We 
assume L 2 V(4>) ~ —12 + L2 ^ f2 0 2 for 0 —> 0 O = 0 which results in A(0 O ) = whereby 

we have set LA = 1 [13] and, as usual, L 2 M 2 = A(A — 4). We consider 2 < A < 4. 
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From s and T, the pressure follows as 

,, , /'*" 7 dT(i„) - 

PWh) = / d(j>H —~- s{(/)H), (7) 

J oo dcpH 

where p(oo) = 0 holds if V'/V < 2 a/2/3 with V'/Vl^oo —> const, corresponding to a 
’’good” singularity at 0 = oo [Hj. We consider only such cases. 

Besides of a proper adjustment of the dilaton potential V (0) to the equation of state, the 
model parameters G 5 /T 3 and L must be fitted, too. Since a direct mapping procedure of 
an input equation of state to the potential is not at our disposal, we use as trial ansatz 


, . V' D + Si0 3 for 0 < 0. 

Ud(0) = - = < 

[ (f, tanh(t 2 0 - t 3 ) + t 4 ) (l - C osh ( ^-fe 3 )0 for 

(demanding differentiability of vd at 0 m hxes L 2 M 2 and Si) and find 


fit to 

0m 

t\ 

t'2 

t’3 

U 

h 

b 2 

h 

g 5 /l 3 

v 2 . 

0.2163 

0.6453 

0.4988 

0.0845 

0.0286 

0.4842 

2.5020 

3.9887 

0.4544 

s/T 3 

0.2430 

0.6480 

0.5023 

0.0855 

0.0344 

0.4844 

2.6162 

4.1458 

0.4586 


with LT C = 1.8036 (fit to v 2 from [IT]) or LT C = 0.5051 (fit to s/T 3 from [Hj). This ansatz 
obeys the Chamblin-Reall IR behavior L 2 R(0 —> oo) ~ e ( 4 i+*4)0_ xhe approach belongs to a 
similar class of holographic models as the model class B in [15] : it has no confinement in the 
sense of pH) for t\ + £4 < \j2j3 and no explicit fermionic degrees of freedom. Our ansatz 
is meant purely to match lattice QCD thermodynamics data in a restricted temperature 
interval. It is therefore interesting to see in future investigations, e.g. whether the explicit 
account of quarks has a similar impact on 0/s as found in the present setting. To set a scale, 
we determine T c in the holographic model by the inflection point of s/T 3 as a function of 
T, and T c = 155 MeV |[TT] is used in the lattice QCD data [H]. The resulting velocity of 
sound squared, v 2 = 1 , the scaled entropy density, s/T 3 , the scaled 

pressure, p/T 4 , and the scaled interaction measure I/T l = (sT — 4 p)/T { are exhibited in 
Fig. 1 together with the lattice QCD data mm- The solid blue (dotted red) curves are 
our optimum fits of v 2 (s/T 3 ) with the parameters of if Circles depict the respective 
quantities at T c . One observes that the softest point, i.e. the minimum of v 2 as a function 
of T/T c , is slightly below unity (see upper left panel in Fig. 1) and the maximum of the 
interaction measure is a little bit up-shifted in comparison with the lattice QCD results (see 
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lower right panel in Fig. 1). One observes also some other minor imperfections of our fits, 
in particular at the lowest temperatures covered by the lattice QCD data, and also at large 
temperatures for the interaction measure. 




FIG. 1: Velocity of sound squared v 2 (upper left panel), scaled entropy density s/T 3 (upper right 
panel), scaled pressure p/T 4 (lower left panel), and scaled interaction measure I/T 4 (lower right 
panel) as functions of T/T c . Solid blue curves: fit to v 2 data, dashed red curves: fit to s/T 3 data, 
circles: position of thermodynamic quantities at T c . Lattice QCD data: black plusses from HU, 
green crosses from [12] : for both sets we use the pseudo-critical temperature T c = 155 MeV. 


In the present setting, the ratio of shear viscosity to entropy density is rj/s = 1 /( 47 t) [15] . as 
usual for the Hilbert action on the gravity side [13 HE]. The ratio of bulk to shear viscosity 
can be calculated via the Eling-Oz formula [18] 


C = / dlog g \ ~ 2 = / 1 dlogT \-2 
r 7 cj, H V d(f>H ' v v 2 dcj)H ' 


( 10 ) 
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or, equivalently |T9] , via the Gubser-Pufu-Rocha formula m 

C v'(M 2 i 

V <t>H V(4 >h) 2 |hn(0o)| 2 ’ 

where hn(4> o) is extracted from solving the perturbation equation 


( 11 ) 


XV 


^ + + » (2 + - 4fl + r - fl + = o, ( 12 ) 


4X V 


Y 


X 2 


3 V' 


4X V 


with initial conditions hn((j)H — e) = 1 and h' u ((j)H — e) = 0 for e —>• 0. The result is 



FIG. 2: Scaled bulk viscosity, (/T 3 , (left) and ratio of bulk to shear viscosity, (/r], (right) as a 
function of T/T c . Line codes as in Fig. 1. 


exhibited in Fig. 2. Remarkable is the reduction of (/rj by 50% at T c in comparison with 
the SU(3) gluon plasma (YM) considered in |20j . To understand this difference, recall the 
adiabatic approximation of [fl3] : X (0) w — • In this approximation, the non-local term 

in becomes unity, h u = 1 (since the coefficient of h u in ( [l2| vanishes, see also [21]), 
and, comparing the values of V' jV ~ 0.6 (0.8) for the QGP (pure glue, see below) at T c 
we hnd the ratio ((C/r]) QGP /(C/rj) y AI )(T C ) ~ 56% (cf. also [15], .22] for recent holographic 
calculations of transport coefficients). 

On the other hand, for (/T 3 , the situation is reversed: at T c , the QGP value is 50% larger 
in comparison to the gluon plasma case; the peak of (/T 3 is located at a larger value 
T/T c fa 1.3. This difference between (/r / and (/T 3 for QGP and for the gluon plasma can 
be attributed to the different number of degrees of freedom as reflected by the scaled entropy 
density s/T 3 . 
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III. CROSS-OVER VS. FIRST-ORDER PHASE TRANSITION 


Remarkably, the ansatz (J8]) for the dilaton potential is the same as used in J2j0j] to describe 
the SU(3) gluon plasma, which displays a first-order phase transition. (This is actually not 
so surprising, as [13] has demonstrated that a two-parameter ansatz for the potential allows 
either for a cross-over, or a first-order phase transition or a second-order transition, depend¬ 
ing on the choice of the parameters. Also, [15] uses a unique ansatz with two parameter sets 
to arrive at a first-order phase transition or a cross-over.) To elucidate the origin of such a 
difference we exhibit in Fig. 3 a few relevant quantities of both optimized models. 




FIG. 3: V'/V (solid black curves), v 2 (dashed red curves, dot-dashed curves are for the adiabatic 
approximation) and 0.1T/T C (blue dotted curves) as functions of 4>h • Left panel: for the pure 
gluon plasma (the grey band covers the unstable and metastable regions), right panel: for QCD 
quark-gluon plasma (fit to v 2 ). Vertical dotted lines bracket the fit range to the lattice data. 


For the pure gluon plasma (left panel) the quantity V'/V has a first maximum of about 0.8 
at 4> h ~ 5.1. In the adiabatic approximation |T3] the velocity of sound squared is 


i.e. a local maximum (minimum) of V'/V is related to a local minimum (maximum) of v 2 . 
If V'/V is sufficiently large, v 2 can go to zero. In fact, the adiabatic approximation is quite 
accurate (compare the red dot-dashed curve (adiabatic approximation) and the dashed curve 
(exact result) in Fig. 3 - left panel). That implies, lifting V'/V sufficiently causes a first-order 
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phase transition, here signalled by v 2 = 0. The entropy density s((pn) is a monotonously 
dropping function, as holds true for the considered examples (and is assumed to hold in 
general in the thermodynamically stable phase, see p3])- Hence, v 2 = 0 corresponds to an 
extremum of T((pn), which is a minimum (maximum) if dv 2 /d<pH < 0 (dv 2 /d(pn > 0). Thus, 
if V'/V is adjusted such that v 2 ((Ph) becomes negative in some (pH interval and then, for 
larger (pH, rises to become positive again, the local minimum of T((pn) is followed by a local 
maximum. (These extrema are very shallow in the left panel of Fig. 3 and hardly visible on 
the used scale.) Such a behavior of T((p H ) leads in turn to the usual loop structure of p(T), 
characteristic of a first-order phase transition. For the case at hand, the pressure is always 
positive. (In contrast, the IHQCD model p3j has one global minimum of T((pn) which gives 
rise to the high-temperature branch and the unstable section of p(0#); the low-temperature 
branch is represented by the line p = 0 corresponding to the thermal gas.) 

Inspection of the same quantities for our fit of the QCD equation of state (see right panel in 
Fig. 3) reveals V'/V < 0.8 everywhere in the considered range of (p from 0 (UV region) up 
to 10 (towards the IR region), that is v 2 > 0 everywhere. Also here, the adiabatic approxi¬ 
mation is quite accurate (compare the red dot-dashed curve (adiabatic approximation) and 
the dashed curve (exact result) in Fig. 3 - right panel). 

Note that Fig. 3 also unravels some uncomfortable features of the ansatz (J8]) with parameters 
adjusted to the lattice Yang-Mills equation of state as in [20] : To catch the shape of thermo¬ 
dynamic quantities in the temperature range (0.7 — 10)T C , the ansatz (J8]) does not qualify to 
continue towards the deep IR region, since, e.g., v 2 becomes negative for <p H > 9.5, signalling 
the break-down of the ansatz’s capabilities. (From the IHQCD viewpoint such a behavior 
is admissible: the point where v 2 = 0 would signal a Hawking-Page phase transition to the 
p = 0 phase, and desirable: the model becomes zero-T confining [T4]. In contrast, our ansatz 
Q is an ad hoc construction to mimic the Yang-Mills equation of state for T > 0.7T C (up to 
10T C ), corresponding to (pH < 8.5 (down to <pn ~ 1-5). It can be supplemented by further 
terms becoming relevant for (pH > 8.5. Thus, it is meaningless to derive from (J8| properties 
of the boundary theory in the IR region.) In contrast, for the QCD parameter adjustment 
(see (|9|) ), the ansatz (J8| seems to be applicable towards the deep IR region. 

Upon an integration of V'/V the potentials V((p) emerge, displayed in Fig. 4. In contrast 
to V'/V, the potentials look quite featureless, both in the region where the softest point 
(minimum of v 2 ) appears for the QCD equation of state (depicted by the arrows) and in the 


region of the first-order phase transition for the Yang-Mills equation of state (grey polygon). 





FIG. 4: The dilaton potentials V (d) in units of — 12L 2 for the ansatz Q with parameter sets from 
© (solid blue or red dashed curve for the fit of u 2 or s/T 3 ; the arrows point to the location where 
Vg has the minimum). The potential ^ with parameters adjusted to the Yang-Mills equation of 
state [20] is exhibited by the dotted black curve (the un/metastable region is depicted by the grey 
polygon). Vertical dashed (dotted) lines bracket the fit range to QCD [TYJ (Yang-Mills [ID]) lattice 
data. 


IV. DISCUSSION AND SUMMARY 

In contrast to the IHQCD model, which covers quite a lot of QCD features both for the pure 
gluon plasma [T|| [23] and for QCD in the Veneziano limit pM], at finite as well as at zero 
temperature together with a direct account of the two-loop t’Hooft running coupling, we 
consider here a simple holographic gravity-single-dilaton model without any explicit a priori 
scale setting. All parameters are adjusted to finite-temperature lattice QCD thermodynam¬ 
ics in a selected temperature range. We formulate a simple criterion to see already at the 
dilaton potential (actually its scaled derivative) whether a first-order phase transition can 
emerge, as for the pure gluon plasma, or a cross-over is encoded, as for QCD at /i = 0. While 
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our focus is clearly on features in a limited temperature range at T > T c , also some section 
of the low-temperature region can be successfully accomodated in the model, leaving the 
deep IR region for further studies. We also stress that we do not require a specified behavior 
of the model outcome in the UV region. Note here that the influence of both asymptotic 
regimes on the equation of state in the considered temperature interval is fairly small: as 
shown in [20] the influence of the UV region on dimensionless thermodynamic quantities 
should not exceed a few percent; for T > 0.7 T c , the deep IR region contributes to p and 
/ as a small integration constant, while s, T, v^, and the viscosities are independent of it. 
Hence, although our dilaton potentials ignore QCD features at T —> 0 and T —> oo, we 
argue that they qualify for further investigations. For instance, supplemented by a fit of the 
quark number susceptibility one can repeat the analysis of [8] with an up-to-date input to a 
holographic study of the phase diagram. Even prior to that we note the interesting drop of 
the ratio (/p by 50% at T c when including quarks. 

Another obvious extension of our studies would be the inclusion of a field dual to the chiral 
condensate (qq), which is responsible not only for the breaking of conformal invariance 
in addition to the gluon condensate as expressed by the trace anomaly but, even more 
importantly, it is in the chiral limit an order parameter of chiral symmetry breaking in 
QCD. Extensive investigations in this direction, albeit for the Veneziano limit QCD, were 
performed in [23], where chiral symmetry breaking is realized by tachyon dynamics. 

In summary, we present an adjustment of a single-field dilaton potential to recently confirmed 
lattice QCD thermodynamics data in the temperature range (0.7 — 3.5)T C . A criterion is 
delivered for ensuring a cross-over at the softest point. 
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